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Abstract 

A new method for characterizing the deformable porous materials with 
non-critical adsorption probes is proposed. The mechanism is based on a 
driving the adsorbate through a sequence of constrained equilibrium states 
with the insertion isotherms forming a pseudo-critical point or a van der 
Waals-type loop. In the framework of a perturbation theory and Monte 
Carlo simulations we have found a link between the loop parameters and the 
host morphology. This allows one to characterize porous matrices through 
analyzing a shift of the pseudo-critical point and a shape of the pseudo- 
spinodals. 

I. INTRODUCTION 

One of the methods for characterization of porous materials^'^ is a retrieving of pore 
size distributions through adsorption probes. Capillary condensation in macro- or meso- 
pores permits to detect them rather reliably because the condensation pressure and tem- 
perature are very sensitive to the pore size. Nevertheless, it becomes more and more 
questionable^ that the true equilibrium behavior is observed in experiments. One of the 
major problems is a well developed metastability leading to very long equilibration times. 
Therefore, the condensation is usually detected as an onset of hysteresis. Nevertheless, 
it has recently been argued"^ that the hysteretic behavior is not necessarily associated 



with an underlying equilibrium phase transition. In the case of micro-pores the situation 
is even more difficult because such pores are too small to support any sharp (bulk- like) 
phase transition. This makes one to search for alternative characterization mechanisms 
which are not based on the critical behavior. In the case of deformable matrices^, such as 
aerogels^'^ or amorphous intercalation compounds^, the problem is additionally compli- 
cated by variations of the host morphology with the guest density or pressure. Thus the 
guest response (e.g., adsorption isotherms) is a superposition of at least two contributions 
- the adsorbate thermodynamics under the host confinement and the adsorbate-induced 
host deformation. 

In this paper we propose a method for a pore size determination through adsorption 
of non-critical fiuids. The main idea is to find an explicit relation between the density- 
dependent matrix morphology and some peculiarity in the fluid thermodynamic response. 
The mechanism is based on our finding^ of negative compressibility states in confined 
dilatating systems. These states can appear as a consequence of a driving path - ei- 
ther a sequence of controlled lateral dilatations^ or a sequential adsorption of prescribed 
portions^^. This imposes an additional constraints^ on the fiuid equilibrium. As a result, 
at certain conditions one can observe a sequence of states which form a van der Waals-type 
loop in the isotherm. Then the loop parameters can be taken as the required specificities 
for a detection of the confining geometry. 

II. MODEL 

The matrix is considered as a collection of slit-like pores with different widths hi 
distributed according to some probability density. For simplicity the pores are assumed 
to be uncorrelated. This allows us to focus on a single pore of fiuctuating width h 
and surface area S. The width is known only statistically, that is, from a probability 
distribution f{h). The sht itself is a part of the hosting system, whose couphng to the 
guest species changes the sht geometry (see below). The fiuid is modeled as hard 
spheres of diameter a = 1, injected into the pore. The hard core potential allows us to 



mimic the excluded volume effects and ensures that the fluid is off-critical (we do not 
consider freezing-type transitions). 
For any h the total Hamiltonian is 

H = Hff + Hfw, (1) 

where Hff is the hard sphere Hamiltonian, Hf^j is the slit potential 



N 
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l/2<Zi<h-l/2. (2) 



We do not take into account a short-ranged fluid-wall attraction, responsible for the 
surface adsorption or layering effects. The inverse-power shape for Hf^ is chosen as a 
generic form of a soft repulsion, with k controlling the softness. In some sense this mimics 
the hydrophobic or non-wetting effects that, as we will see later, are essential. On the 
other hand, this simple model takes into account that the walls are impenetrable and 
ignores possible localized adsorption effects, which could complicate the overall picture. 
For technical purposes we are working with A; = 3. Moreover, our results are quahtatively 
insensitive to a particular choice of k. 

In principle, one may consider two adsorption mechanisms. One is an equilibrium 
between the pore and the bulk fluids. In that case would change according to the 
chemical potential /j, = /lb. In this paper we treat another mechanism. Namely, an 
adsorption in doses^*^, when N is injected in prescribed portions. In that case we have 
a sequence of the fluid equilibrium states corresponding to different N. The sequence 
is determined by an injection procedure, which imposes an additional constraint. It 
should be emphasized that choosing the model we tried to retain only a minimal number 
of relevant ingredients: the absence of liquid-gas transition in the fluid bulk, excluded 
volume effects, and a density-dependent heterogeneous confinement. 

III. PERTURBATION THEORY 
A. Conditional insertion isotherm 



At any pore width h the free energy can be represented as 



(3) 



where Fo{h) is the free energy of a reference system, and (...)o is the average over the 
reference state. Taking a spatially confined hard sphere system (the one with ^4 = 0) as 
a reference, we consider a first order perturbation^ for the conditional free energy 



/3F{h) = /3Fo{h) + /3Sp{h) I n dziHf, 



(4) 



where /? — l/{kT) and p{h) is the pore density in the "slab" approximation 

p{h) ^ 



(5) 
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ignoring a non-monotonic behavior of the density profile with increasing pore density. 
This is reasonable for wide pores and low fluid densities. The reference part is estimated 
in the excluded volume approximation, while the perturbation contribution is N'^{h), 
such that the total conditional free energy is 

"1 - hp{h) 



j3F{h) = -TV In 
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where h = 27r(T^/3 is the excluded volume factor, is the thermal de Broglie length and 



■^{h) = IQA- 
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The conditional insertion isotherm can be calculated in the standard way 
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which relates the fluid density p{h) and the chemical potential p{h). 



B. Matrix dilatation 

As is discussed above, the pore dilatation can be taken into account, assuming that 
the width h is density-dependent. Real insertion materials are usually rather complex 
(multicomponent and heterogeneous). For that reason one usually deals with a distribu- 
tion of pore sizes or with an average size. In this context we assume that h is known only 



statistically and the probability distribution f{h\p) is conditional^ to the guest density p, 
which should be selfconsistently found from 

p = J dhf{h\p)p{h). (9) 

This relation reflects an injection sequence p{h) and the matrix dilatation. Note that 
the average density p is not uniquely determined by the particles insertion but is also 
conditional to the matrix response. Namely this constraint makes the problem non-trivial, 
inducing a coupling between the insertion and the dilatation modes. Taking our result 
for p{h), we can focus on the insertion isotherm, averaged over the width fluctuations. 

p{p) = J dhf{h\p)p{h). (10) 

This, however, requires a knowledge on the distribution f{h\p). Even without resorting 
to a concrete form for f{h\p), it is clear that the matrix reaction can be manifested as 
a change in the distribution width or/and the mean value. One of the simplest forms 
reflecting at least one of these features is a 5-like distribution, ignoring a non-zero width: 

f{h\p) = 6[h-h{p)], (11) 

where S{x) is the Dirac 5-function and the mean pore width h{p) is density-dependent. 
For concreteness we consider the swelling behavior^. Therefore, the mean pore width 
h{p) increases with p. For illustration purposes we take the (S'-shaped form, considered 
previously^'^ 

h{p) = ho{l + tanh[A(p - po)] + tanh[Apo]). (12) 

This form resembles a non-Vegard behavior, typical for layered intercalation compounds^. 
The dilatation is weak at low densities (p << po), the most intensive response is at 
p ^ Po, and then the pore reaches a saturation, corresponding to its mechanical stability 
limit. Here A is the matrix response constant or dilatation rate, controlling the slope 
near p ^ pg, and ho is the mean pore width in the absence of insertion. This should be 
considered as a generic form that samples a non-linear increase of the average pore width. 



resolved at a given density p. The model mimics certain heterogeneity aspects which, as 
we believe, are essential. On the other hand, this choice allows us to map the problem in 
question onto the negative compressibility problem^ in systems with fluctuating geometry. 
Prom eq. (9) the average density is found to be 

Changing the surface density N/ S (by increments in N) we vary the average pore density 
p. This allows us to ehminate N/S in the favor of p in all thermodynamic functions. 
Combining eqs (8), (9) and (10) we obtain 



pp{p) = In 
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In the case of wide pores and a weak dilatation we expand j3p in. terms of I/Hq and A, 
obtaining a generic van der Waals form 



I3p{p) = In 



A3p 
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Therefore, as in the previous study^, we have an interplay of several effects - the packing 
(first two terms), the fluid- matrix interaction (density-independent), and the dilatation 
(linear in density) . It is seen that a coupling of the fluid-slit repulsion (A) , spatial confine- 
ment (ho) and the dilatation (A) acts as an effective infinite-range fluid-fluid attraction. 
The latter induces an inflection point (or even a loop) in the isotherm p{p). Introducing 
a dimensionless temperature T* = 1/(4^4*), and solving 

we find the pseudo-critical parameters at which the inflection point appears. 

1 * 4 A 

^^=36' ^^=276^ ^^^^ 

Plugging these back into eq. (14) we get the third parameter pc = l^'cipi ho, A). 

As has already been mentioned, the loop is not a signature of a phase transition. It is 
a sequence of equilibrium states, constrained by a coupling between the insertion and the 



dilatation modes. The physical mechanism is quite similar to the previously reported^. 
We deal with a competition of two effects: the interparticle and wall-particle repulsions 
tend to increase the insertion energetic cost /i with increasing p, while the pore dilatation 
diminishes it. As a result, at certain parameters the isotherm becomes nonmonotonic. For 
a purely attractive pore {A < 0) the inflection point does not appear. Therefore, the wall- 
particle repulsion is an essential ingredient. This suggests that a suitable combination of 
the repulsive and attractive parts is a criterium for choosing the adsorption probe. As 
it should be, the inflection (loop) feature disappears in the bulk limit /iq — oo or in the 
case of rigid matrices A ^ 0. 

In order to improve the accuracy of the theoretical results, the reference part is replaced 
by the Carnahan-Starling form 



This allows us to avoid the unphysical behavior at high densities and escape from the 
flaws of the excluded volume approximation (14). In particular, all the pseudo-critical 
parameters {pc, T*, pc) appear to be related to the parameters of the pore width distribu- 
tion {hQ,pQ,A). Therefore, driving the system close to the pseudo-critical point permits 
to estimate the density dependent pore size distribution. 



In order to verify the existence of the loop predicted by the perturbation theory as 
well as to test the accuracy of the theoretical estimation, the simulation of the model in 
question has been carried out appling the Widom test particle insertion method^^. As 
in our previous work^ we applied the conventional canonical NVT MC simulations of a 
confined hard sphere fluid. The simulation cell was parallelepiped in shape, with parallel 
walls at surface separation h, and constant surface area S — x Ly — 12 x 12. The 
periodic boundary conditions were applied to the X and Y directions of the simulation 
box; the box length in the Z direction is flxed by the pore width. For a given pore 





(18) 



IV. SIMULATION 



width the adsorbed fluid density is chosen according to the Eq.(5) varying the number 
of particles, N. At frequent intervals during this simulation, a coordinate r^^i has been 
generated randomly and uniformly over the simulation cell. For this value of r^v+i, we 
have computed (exp(— /3if(rjv+i)))Ar averaged over all generated trial positions. In such 
a manner the chemical potential of the system with density p has been calculated. 
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_(exp(-/3i7(riv+i)))iv. 
where H is the system Hamiltonian. 

Each simulation runs 4 x 10^ MC cycles, whith the first half for the system to reach 
equilibrium whereas the second half for evaluating the ensemble averages. 



V. RESULTS 

Adsorption isotherms are plotted in Figure 1. It is seen that the infiection point around 
p ~ 0.3 transforms into the loop with increasing repulsion A (or decreasing temperature 
T*). As expected^, the perturbation theory underestimates the magnitude of A at which 
these effects appear. The loop gets sharper with increasing dilatation rate A and becomes 
more pronounced with decreasing pore width ho. It is remarkable that the fitting on the 
pseudo-critical (infiection) point (part (a)) or on the loop (part (b)) allows us to find the 
set of the distribution parameters {ho, A, po) in a perfect agreement with the simulation. 
Attempts to fit under other criteria (e.g. the least-square deviations) give incorrect values 
of these parameters. 

In agreement with our theoretical prediction (17), the simulation results confirm that 
the pseudo-critical point (and the loop) is reachable at any appropriate choice of the model 
parameters. This is in contrast to what was reported^ for the negative compressibihty 
states in the same model with lateral stretching (increasing S while N is fixed) . There the 
loop has been found to appear only if the sht reaction A reaches some threshold value A* 
which involves a combination of ho, po and A. In other words, the variation of the average 
density p due to the transversal dilatation should dominate its variation, induced by the 



changes in the surface density N/S due to the lateral stretch. This difference suggests 
that the driving path (namely, changing the fluid density by a lateral stretch or by a 
particle injection) is really important for this sort of systems. On the other hand, in the 
context of the present study, the fact that the pseudo-critical parameters do not involve 
any additional condition, is quite favorable for making a link between these parameters 
and the matrix morphology. 

In Figure 2 the pseudo-critical parameters are analyzed. It is seen that Pc grows with 
the dilatation rate A almost linearly exhibiting a characteristic kink. The latter signals 
a crossover from a weak to a strong dilatation regime. The kink shifts to lower A and 
becomes less pronounced with decreasing ho (that is why it is almost undetectable for 
ho = 5). At low A the pseudo-critical density is practically independent of ho, while at 
high A it approaches po- In agreement with our simplified estimation (17), T'^*(A) can 
be reasonably approximated by a straight line (the lower inset). Nevertheless, it exhibits 
a weak inflection point at the same A (around A = 4) as pc does. The pseudo-critical 
temperature incresases with decreasing ho- This is also coherent with the estimation (17) 
predicting a decrease in the slope inversly proportional to ho- In contrast, the insertion 
energy pc decreases and becomes insensitive to ho with increasing A (the upper inset). 
This agrees well with the physical intuition that the pore dilatation diminishes an overlap 
between the repulisive wall flelds, making the insertion less energetically consuming. 

In practice, an exact capture of the inflection point could be too demanding. Driving 
at rather low temperatures (T < T*), such that the system develops the loop, should 
be much easier. In that case a shape of the pseudo-spinodals could also be used as an 
indicator of the matrix morphology. These curves are plotted in Figure 3. It seen that the 
pseudo-spinodals become broader with decreasing dilatation rate and shift to higher p with 
increasing A. Comparing this with the previous figure, we can see again that pc converges 
to Po with increasing dilatation rate A. If necessary, other quantitative characteristics, 
such as the curve asymmetry with respect to Pc or p values at the spinodal points, can 
also be easily extracted. 

These results illustrate a correlation between the pseudo-critical point (or the loop) 



parameters and the density-dependent pore geometry. We did not try to obtain these 
parameters from the simulation because the procedure is time consuming, while new 
physical insights, compared to the anlytical estimations, are difficult to expect. The 
only promising point is to consider narrow {1 < ho < 2) pores, for which the first-order 
perturbative scheme could fail even quahtatively. This point is left for a future analysis. 

VI. DISCUSSION AND CONCLUSION 

A method for characterizing the deformable porous materials by adsorption of non- 
critical fluids is proposed. It is based on a driving the host-guest system through a 
sequence of constrained equilibrium states, forming a loop in the insertion isotherms. 
The constraint appears from a coupling of the insertion and the matrix dilatation modes, 
inducing a competition between the repulsive interactions and their relaxation by the 
host dilatation. As a result the insertion isotherms develop pseudo-critical points or 
loops, whose parameters are shown to be correlated with the parameters of the pore 
size distribution. In particular, we have found that the magnitudes of the geometrical 
parameters {ho, po, A) which reproduce well the inflection point or the loop, perfectly 
agree with their values in the simulation. This suggests that the peculiar points (with 
dji/dp = or dp/ dp = 0) are extremely sensitive to the matrix morphology, while 
the overall isotherm shape is less important (unless on can achieve an exact isotherm 
coincidence). This is coherent with the recent finding^^ that a matching of the inflection 
points (with dp/ dp = 0) correctly reproduces the remaining part of the experimental 
isotherm. From the information-theoretical point of view^^ such a criterium (the inflection 
matching) corresponds to a maximal amount of information one can get from a fitting 
of a model to a given experimental or simulation result. In our case we deal with the 
pseudo-critical point {dp/ dp = 0). Fitting on this peculiarity corresponds to a singular 
information rate^^, that allows us to update the geometrical parameters from their initial 
estimation in a limited range of p (around the singularity). Namely this type of singularity 
is used to characterize porous media^ through a shift of the fluid condensation or freezing 



point. 

Although we have taken a simphfied S-like distribution, our conclusion should hold in 
general because the characteristic loop feature is directly related to the dilatation mode. 
The latter translates into a shift of the mean pore size, that has been taken into account. 
Nevertheless, the role of a non-zero distribution width remains to be investigated. The 
only restriction is that the number of the detected loop parameters should be (at least) 
equal to the number of parameters determining the pore size distribution. This restriction 
is not too strong because in practice the distribution is usually characterized by only few 
parameters (e.g. leading mean pore sizes and dipersions around them). On the other 
hand, the proposed method is mainly intended for a complementary detection of the 
porosity range (e.g. micropores) where other techniques fail. For illustration and technical 
purposes we were working with rather large pores 10 < h{p) < 30. This, however, does 
not affect our main concept. 

It should be emphasized that, since we deal with a constrained equilibrium, the loop 
does not violate the thermodynamic stability^^. Therefore, this feature can be readily 
found under suitable experimental conditions, similar to those realized in the low-pressure 
adsorption into microporous matrices^^. In this respect it is worth noting that quite similar 
loop features have been detected^" for nitrogen and argon adsorption in microporous 
carbons and synthetic zeolites. As in our case, the fluids have been inserted in sequential 
doses. Nevertheless, in contrast to our assumption, there is no evidence for an insertion- 
induced dilatation of these slit-shaped materials. The loops have been interpreted as a 
signature of a monolayer- induced micropore filling^^. In this context we can propose an 
alternative interpretation of our results, making a link to the experimental observations. 
Namely, it is well known^'^ that the micropore filling starts from the narrowest pores 
{ho in our laguage). As the fluid pressure (or density) increases, the larger micropores 
progressively fill. Microscopically this effect appears from an interplay of repulsive and 
attractive parts of the wall-particle interaction. The latter depends on the chemical 
nature of the probe and the matrix geometry and composition. In order to extract general 
features, one may accept this effect as an experimental fact, trying to see its consequences. 



Therefore, in a given density domain the fluid resolves a quite narrow range of pore 
sizes because the narrow pores are aheady blocked by a preceding dose, but the flUing 
of larger pores is not thermodynamically favorable at this density. This fact is captured 
by our distribution f{h\p) that can be re-interpreted as an effective density-dependend 
matrix heterogeneity, with A being the heterogeneity index. In other words, an increase 
of the mean pore size h{p) may result from a pore blocking effect. This interpretation is 
coherent with the one proposed by the authors of^°, who claimed that the S-shaped feature 
(the loop) signals a switching from the ultramicropore to micropore flUing. Nevertheless, 
this qualitative analogy between the matrix swelling and heterogeneity requires further 
analysis, including a competition of the repulsive and attractive interactions. This issue 
is left for a future study. 
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FIGURES 

FIG. 1. Insertion isotherms for a weakly (a) and strongly (b) repulsive slit. The other 
parameters are ho = 10, A = 15, po = 0.3 (both theory and simulation) 

FIG. 2. Variation of the pseudo-critical parameters with the dilatation rate A. po = 0.3, 
ho = 5-squares, ho = 10-circles 

FIG. 3. Pseudo-spinodal curves at different dilatation rates A. The other parameters are 
^0 = 7, Po = 0.3 
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